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Abstract 

The dynamics of two 1/2-spin qubits under the influence of a quantum Heisenberg XY type 
spin-bath is studied. After the Holstein-Primakoff transformation, a novel numerical polynomial 
scheme is used to give the time-evolution calculation of the center qubits initially prepared in a 
product state or a Bell state. Then the concurrence of the two qubits, the z-component moment of 
either of the subsystem spins and the fidelity of the subsystem are shown, which exhibit sensitive 
dependence on the anisotropic parameter, the temperature, the coupling strength and the initial 
state. It is found that (i) the larger the anisotropic parameter 7, the bigger the probability of 
maintaining the initial state of the two qubits; (ii) with increasing temperature T, the bath plays 
a more strong destroy effect on the dynamics of the subsystem, so does the interaction g$ between 
the subsystem and the bath; (iii) the time evolution of the subsystem is dependent on the initial 
state. The revival of the concurrence does not always means the restore of the state. Further, the 
dynamical properties of the subsystem should be judged by the combination of concurrence and 
fidelity. 
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I. INTRODUCTION 

Solid-state devices, in particular, ultra-small quantum dots []| with spin degrees of 
freedom embedded in nanostructured materials, compared with other physical systems such 
as ions in trap {2} , are more easily scaled up to large registers and they can be manipulated 
by energy bias and tunneling potentials 3J. The key building block of a quantum processor 
consists of two entangled quantum bits. Thus the spin system is one of most promisingcan- 
didates for quantum computation owing to long relaxation and decoherence times 



However, the spin qubits are open systems which is impossible to avoid interactions with 
their environments flflflfl. Fpifc the states of the q ubits will reha* into a set of "pointer 
states" in the Hilbert space [lOj; and the entanglement between the spin qubits will also 



vanish. Yet the entanglement is the most intriguing feature of quantum composite system 



and the vital resource for quantum computation and quantum communication [5|, llll]- These 
are so-called decoherence and disentanglement processes. These two disadvantages will not 
be overcome until the modelling of the surrounding environment or bath of the spin systems. 

For solid state spin nano-devices, the quantum noise mainly arises from the contribution 
of nuclear spins, which could be regarded as a spin environment. Recently, there are lots of 
works were devoted to study the behavior of center spins under the strong non-Markovian 



influence of a spin-bath [li 



131 ] . Lucamarini and co-workers made use of perturbative 
expansion method [ijj and mean-field approximation [ijj to study the temporal evolution 
of entanglement pertaining to qubits interacting with a thermal bath. They found entangled 
states with an exponential decay of the quantum correlation at finite temperature. Hutton 
and Bose 16| investigated a star network of spins at zero temperature, in which all spins 
interact exclusively and continuously with a central spin through Heisenberg XX couplings 



of equal strength. Their work was advanced by Hamdouni and co-workers [171 ]. who 
derived the exact reduced dynamics of a central two-qubit subsystem in the same bath 
configuration. And they also studied the entanglement evolution of the central system. 
Yuan and co-workers j^] used a novel operator technique to obtain the dynamics of the two 
coupled spins in quantum Heisenberg XY high symmetry spin model. The results of all the 
above works are very exciting. Yet their methods are of some kinds of complex analytical 
derivations. And in Ref. [21] , their analytical results are dependent on some particular initial 



2 



states and essentially the interaction between the spins in their model is isotropic. Here we 
introduce a "half analytical and half numerical" method to solve such an open quantum 
system problem in an anisotropic Heisenberg XY model. The present model involves the 
Heisenberg XY interaction that has broad applications for various quantum information 



processing systems, such as quantum dots, Cavity-QED, etc[18|, [19|, |20|, |21[. Besides, our 
method is initial states independent. 

In this paper, we study an open two-spin-qubit system in a spin bath of star-like config- 
uration, which is similar to the cases studied in Ref. 0, \v\. But the two qubits' distance 
is far enough so that the direct coupling between them could be neglected. Then we can 
concentrate on discussing the role of the bath in this model. First, we use Holstein-Primakoff 
transformation to reduce the model to a effective "spin-boson" Hamiltonian. Then we ap- 
ply a numerical simulation to obtain the reduced dynamics of the two-spin qubits. During 
our numerical calculation, there are no approximations assumed and the initial state of the 
subsystem (consisted by the two spin qubits) can be arbitrary. It is well-known that the 
concurrence is a measure of entanglement degree between two spin qubits and the fidelity is 
also an important property, that has been widely applied into quantum coding theory 22] . 
Thus some results about these quantities in the thermal limit will be given in the latter part 
of this paper. The rest of this paper is organized as following. In Sec. HI] the model Hamil- 
tonian and the operator transformation procedure is introduced. In Sec. IHIl we explain the 
numerical techniques about the evolution of the reduced matrix for the subsystem. Detailed 
results and discussions can be found in Sec. [TV] We will conclude our study in Sec. [V] 

II. MODEL AND TRANSFORMATION 

Consider a two-spin-qubit subsystem symmetrically interacting with bath spins via a 
Heisenberg XY interaction: both the subsystem and the bath are composed of spin-1/2 
atoms. Every spin in the bath interacts with each of the two center spins of equal strength, 



similar to the cases considered in 



231 ] . The Hamiltonian for the total system is 
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divided as three parts: 



H = H S + H SB + H B . (1) 
H s = /i « + ^), (2) 



N 

9o 



H ™ = VW E K 1 + + °S»°?) + ( X - + , (3) 

2ViV ~T 



i=i 

iV 

^ = 2^ E K 1 + W + ^ - ^M] ■ ^ 





Here, Hs and H B are the Hamiltonians of the subsystem and bath respectively, and Hsb de- 
scribes the interaction between them (^,[^,0]. /xo represents the coupling constant between 
a locally applied external magnetic field in the z direction and the spin qubit subsystem. 
li — 1 < 7 < 1 is the anisotropic parameter. When 7 = 0, it is of an isotropic case [9|. In 
the following part of this paper, we only talk about cases with positive 7 for the symmetry 
of the spin star structure, cx^, crj^ and cr^ (z=l,2) are the operators of the qubit subsystem 
spins, respectively. By Pauli matrix, the operators read 

•1 1 "•=! °~ ! 1. '-=[ : 1 (5) 

of and of are the corresponding operators of the ith atom spin in the bath. The indices i 
of the summation for the spin bath run from 1 to N, where N is the number of the bath 
atoms, go is the coupling constant between the qubit subsystem spins and bath spins, 
whereas g is the coupling between the bath spins. 

Using o x = (a + + o~),o y = —i(a + — o~), we can rewrite Hamiltonians. (j3J) and (0J as: 

N N 

'N 

1 

+ E a t^°02 + O-02) + E a i( a 02 + 70Q2) 
i=l i=l 

N 

H b = jj E + °7°J) + ("t*I + °7°t)\ • ( 7 ) 

*¥=3 

Substituting the collective angular momentum operators J± = J2iLi °t Eqs. (jSJ) 



TT 90 

nsB ■ 



E a ^^ a oi + ff oi) + E a * ( a oi + 70r oi^ 

L i=l i=l 

N N 



(6) 
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and (JTj), we get 

H S B = ~J= [J+(l^01 + ^Ol) + J -( a 01 + 7^0l) + J +(l°02 + °w) + J -(°02 + 7^02)] , ( 8 ) 

H B = {J+J+ + J-J-) + 04 J- + J-J+- 2j)\ . (9) 

where j = N/2. After the Holstein-Primakoff transformation 25j, 

J + = b+(y/2j-b+b), J- = {^2j-b+b)b, (10) 

with [6, 6 + ] = 1 and in the thermodynamic limit (i.e. iV — > 00) at finite temperatures, the 
Hamiltonian, Eqs. jS} and (Q, can finally be written as 

H S b = go [b + (icr+ + a- x + 7( t+ + 0-02) + b(p$ x + 7^01 + °"02 + 7^ra)] > ( n ) 
^ = ^[ 7 (& +2 + 6 2 ) + 26+6]. (12) 

The transformed Hamiltonian describes two qubits interacting with a single-mode thermal 
bosonic bath field, so the analysis of the model is just like a nontrivial problem in the 
field of cavity quantum electrodynamics Q, Q • We note here that due to the transition 
invariance of the bath spins in our model, it is effectively represented by a single collective 
environment pseudo-spin J in Eq. (|T0l) . After the Holstein-Primakoff transformation and 
in the thermodynamic limit, this collective environment pseudo-spin could be considered 
a single-mode bosonic thermal field. The effect of this single-mode environment on the 
dynamics of the two qubits is interesting. In Sec. IIVt we will show some results, for 
example, the revival behavior of the reduced density matrix or entanglement evolution of 
the subsystem spins. This can be used in real quantum information application. 



III. NUMERICAL CALCULATION PROCEDURES 

The initial density matrix of the total system is assumed to be separable, i.e., p(0) = 
® pb- The density matrix of the spin bath satisfies the Boltzmann distribution, that 
is Pb = e~ HB / T /Z , where Z = Tr (e~ Hs / T ) is the partition function, and the Boltzmann 
constant ks has been set to 1 for simplicity. The density matrix p(t) of the whole system 
can formally be derived by 

p(t) = exp(-iHt)p(0)exp(iHt), (13) 
p(0) = Ps (0)®Pb(0), (14) 
p s (0) = |V>(0)>(V(0)|. (15) 
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n order to find the density matrix p(t), we follow the method suggested by Tessieri et al. 



261 ] . The thermal bath state ps(0) can be expanded with the eigenstates of the environment 



Hamiltonian Hb in Eq. (1121) : 



M 



Pb(0) = ^2 l&n^m^ml* ( 16 ) 



m=l 



U m = 7? , (17) 

M 

Z = J2e~ Em/T - (18) 

m=l 

Here \4> m ), m = 1, 2, 3, • • • , M, are the eigenstates of Hb, and E m the corresponding eigenen- 
ergies in increasing order. M is just the number of eigenstates considered in this summation. 
With this expansion, the density matrix p(t) can be written as: 

M 

p(t) = y %2u m \* m (t))(* m (t)\. (19) 

771=1 

Where 

|* m (t)> = exp(-iift)|tt TO (0)) = Z7(t)|* TO (0)). (20) 

The initial state is 

|*m(0)> = |^(O))|0 m ). 

The evolution operator U(t) can be evaluated by different methods. In Ref. [jj, they use a 
unique analytical operator technique. Here, we apply an efficient numerical algorithm based 



on polynomial schemes 27], |28|, |29J into this problem. The method used in this calculation 



is the Laguerre polynomial expansion method we proposed in Ref. [27(], which is pretty well 
suited to many quantum systems, open or closed, and can give accurate result in a much 
smaller computation load. More precisely, the evolution operator U (t) is expanded in terms 
of the Laguerre polynomial of the Hamiltonian as: 

( 1 \ a+l CO y . x k 

L%(H) is one type of Laguerre polynomials {30I as a function of H, where a (—1 < a < 00) 
distinguishes different types of the Laguerre polynomials and k is the order of it. In real 
calculations the expansion has to be cut at some value of /c max , which was optimized to be 
20 in this study (We have to test out a k max for the compromise of the numerical stability 



in the recurrence of the Laguerre polynomial and the speed of calculation). With the 
largest order of the expansion fixed, the time step t is restricted to some value in order 
to get accurate results of the evolution operator. At every time step, the accuracy of the 
results will be confirmed by the test of the numerical stability — whether the trace of 
the density matrix is 1 with error less than 10~ 12 . For longer times the evolution can be 
achieved by more steps. The action of the Laguerre polynomial of Hamiltonian to the 
states is calculated by recurrence relations of the Laguerre polynomial. The efficiency of 
this polynomial scheme 27| is about 9 times as that of the Runge-Kutta algorithm under 
the same accuracy condition used in Ref. [26|. When the states \^ m (t)) are obtained, the 
density matrix can be obtained by performing a summation in Eq. (1191) . 



Although theoretically we should consider every energy state of the single-mode bath 
field: M — > oo, but the contributions of the high energy states \<f> m ),m > m c [rn c is a cutoff 
to the spin bath eigenstates) are found to be neglectable due to their very tiny weight value 
ui m , as long as the temperature is finite. That is to say, the M in Eqs. ( TIB!) to (Tl9l could be 
changed to mc- Then we use the following equation in real calculation: 

mc 

p(t) = Y^u m \V m (t))(V m (t)\. (21) 

m=l 

After obtaining the density matrix of the whole system, the reduced density matrix is 
calculated by a partial trace operation to p(t), which trace out the degrees of freedom of the 
environment: 

p s (t) = Tt b (p(t)) . (22) 

For the model of this paper, ps = is the density matrix of the open subsystem consists 

of two separate spins, which can be expressed as a 4 x 4 matrix in the Hilbert space of the 
subsystem spanned by the orthonormal vectors 1 00) , |01), 1 10) and |11). The most general 
form of an initial pure state of the two-qubit system is 

|^,(0)) = Q |00) + P\U) + 7|01> + 5|10), (23) 
with \a\ 2 + \(3\ 2 + | 7 | 2 + \5\ 2 = 1. (24) 
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IV. NUMERICAL SIMULATION RESULTS AND DISCUSSIONS 



When the reduced density matrix is determined, any physical quantities of the subsystem 
can be readily found out. In the following we will discuss three important physical quantities 
of the subsystem which reflect the decoherence speed, the entanglement degree and the 
fidelity of the subsystem state. These quantities are (i) the moment of spin-01, here we 
choose the first s pin (cqi), which demonstrates the decoherence rate of the system; (ii) the 



concurrence 



31 



32] for the two spins of the open subsystem. The concurrence of the two 



spin-1/2 system is an indicator of their intra entanglement, which is defined as 3JJ: 



C = max{X 1 - A 2 - A 3 - A 4 , 0}, (25) 

where A, are the square roots of the eigenvalues of the product matrix ps(<? y ®cr y )p* s (a y ' ®a y ) 
in decreasing order; (iii) the fidelity [33], which is defined as 



Fd(t) =Tr 5 [pidcai(t)p(t)]. (26) 

Pideai(^) represents the pure-state evolution of the subsystem under H s only, without 
interaction with the environment. The fidelity is a measure for decoherence and depends 
on pideah is equal to one only if the time dependent density matrix p{t) is equal to pideaiif). 
The corresponding results and discussions are divided to two subsections according to 
different initial states. 



For the product states, the present paper focuses on the entanglement generation by the 
spin bath and dose not involve the revival of the initial state for C(t = 0) = 0. Thus in 
the subsection IIV A\ we give out the dynamics of concurrence and Oq X . For the Bell states, 
since the system can evolve to a completely different state from the initial one and has the 
same concurrence C(t) > 0, we should give out the evolution of concurrence and fidelity in 
subsection IIVBI 



A. Product states 



First we show the evolution of concurrence and z component moment as functions of 
anisotropic parameter 7 from two initial product states |11) (to see Fig. [1]) and 1 10} (to see 
Fig. [3]). It is obvious that the entanglement between the two subsystem qubits can arise 
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(a)C(t) 0>)Hx(t)) 

FIG. 1: Time evolution for (a) Concurrence, (b) the moment of spin-01 from an initial two-qubit 
state of |^(0)) = 1 11) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 
(dashed curve), 7 = 0.6 (dot dashed curve), 7 = 1.0 (dotted curve). Other parameters are fio = 2g, 
9o = 9,T = g. 
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(a)C(t) 



(b)(^i(«)) 



FIG. 2: Time evolution for (a) Concurrence, (b) the moment of spin-01 from an initial two-qubit 
state of I -0(0)) = 1 11) at different values of temperature: T = 0.2g (solid curve), T = g (dashed 
curve), T = 5g (dot dashed curve). Other parameters are (jlq = 2g, go = g, 7 = 0.2. 



from the interaction with the bath. Yet this kind of effect of the spin-bath is decreased 
with increasing the anisotropic parameter 7. And the above variation depends on the 



initial states: with \ip(0)) = |11) (Fig. 1(a)), when 7 > 0.87, the concurrence of the two 



qubits will be always kept zero as initialed; with ^(0) = |01) (Fig. 3(a)), the entanglement 
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(a)C(t) (b)(cr z oi(t)) 

FIG. 3: Time evolution for (a) Concurrence, (b) the moment of spin-01 from an initial two-qubit 
state of 1^(0)) = 1 01) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 
(dashed curve), 7 = 0.6 (dot dashed curve), 7 = 1.0 (dotted curve). Other parameters are fio = 2g, 
9o = 9,T = g. 





(a)C(t) 



(b) 



FIG. 4: Time evolution for (a) Concurrence, (b) the moment of spin-01 from an initial two-qubit 
state of IV'(O)) = 1 01) at different values of temperature: T = 0.2g (solid curve), T = g (dashed 
curve), T = 5g (dot dashed curve). Other parameters are [Aq = 2g, go = g, 7 = 0.2. 



can always be created to some extend; and if 7 approaches to (the isotropic case), the 



concurrence can increase as high as 0.8 over some period of oscillations. In Fig. 1(b) and 



Fig. 3(b) with 7 increasing, the oscillation amplitudes of the curves become smaller and 
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gt gt 

(a)C(t) (bX^t)) 

FIG. 5: Time evolution for (a) Concurrence, (b) the moment of spin-01 from an initial two- 
qubit state of |^(0)) = 1 01) at different values of coupling strength between subsystem and bath: 
go = 0.5(7 (solid curve), go = g (dashed curve), go = 2g (dot dashed curve). Other parameters are 
Mo = 2g, 7 = 0.2, T = g. 

smaller, which means that the coherence of the subsystem approaches to be lost. For the 
initial state |11), after the first spin flip (the sign of changes from positive to negative) 
for 7 > 0.6, it can not flip again. However, for \ip(0)) = |10), the spin can flip after a period 
of time even for 7 = 0.6. Therefore, the increase of entanglement depends sensitively on 
the anisotropic parameter. 



The bath is in a thermal equilibrium state at different temperature, which effect is 
shown in Fig. [2] and Fig. HI In these two figures, the anisotropic parameter 7 is kept as 
0.2. We can find that (i) at a very low temperature, OqiCO displays a nearly periodical 
oscillation, which is identical with the two-photon resonance of two two-level atoms in a 
cavity. And the subsystem entanglement can be raised to a comparatively degree; (ii) with 
increasing temperature, the oscillation amplitudes of the curves are damped due to the 
thermal bath. For the concurrence, C(t) — ► means to approach a most separated state 



(to see the dot dashed curve in Fig. 2(a) and Fig. 4(a)). For a^ l7 it means the degener- 



ation of its magnetic moment (to see the dot dashed curve in Fig. 2(b) and Fig. 4(b)). 
Therefore it is clear that the subsystem loses its memory faster as the temperature increases. 
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In Fig. 1(a) and 3(a) , we can find that the entanglement between the two initial separated 
spins can be generated with the help of the single-mode thermal bosonic bath field. Assume 
that the system is initially prepared in |01) (or |10)), a pure state. On one hand, when the 
interaction between the system and spin bath is turned on, one spin could drop from the 
excited state and simultaneously the other spin could jump absorbing the boson just emitted 
by the former. This process induces the entanglement of two spins. On the other hand, it 
also evolves into a mixed state resorting to the bosons provided by the single-mode boson 
field. Thus, the coupling between the system and its environment leads to the entanglement 
between two initial separated qubits. 

Then we keep the bath at a moderate temperature T = lg to find out the effect of the 
coupling strength g in Fig. At a smaller value go = 0.2g, the weak interaction with the 
bath will make both the concurrence and display a pseudo-periodical behavior; on the 
contrary, at a larger value go = 5g, their dynamics is too strongly disturbed by the bath 
to be utilized. Thus in real applications, the coupling between the subsystem and the spin 
bath should be reduced. 



B. Bell states 



In the cases that the subsystem is prepared as a most entangled state C = 1 (Bell 
states), the anisotropic parameter 7 still makes an important effect on the time evolutions 
of concurrence and the Oq X . When ^(0) = l/v^QlO) + |01)), the concurrence (to see Fig. 



6(a) ) of subsystem is always revived to C ~ 0.8 after some time of oscillation at small value 

at 



of 7. The results can be proved by the revival fidelity of the subsystem in Fig. 6(b) 



every summit, the state mainly consists of its initial state. But an interesting phenomena 



is found by the comparison of Fig. 8(a) with Fig. 8(b). It is noticed that two summits 



disappear in the interval of 0.0 < gt < 8.0! Thus we analyze the states at the first three 



summits in Fig. 8(a). It is found that at gt = 2.448, the most component of the subsystem 



state is l/v^GH) — |00)); at gt = 4.960, the state of the two qubits is very near to a 
combination of l/\/2(|ll> - z|00)) and 1/V2(|11) - |00)); at gt = 7.480, the main part of 
the state comes back to its initial state. So the concurrence can not determine the concrete 
state of the subsystem in the present case. Even if the concurrence can be restored, the 
state is not always the same as the initial one. Only the combination of the concurrence 
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gt gt 

(a)C(t) (b)Fd(t) 



FIG. 6: Time evolution for (a) Concurrence, (b) Fidelity from an initial two-qubit state of \tp(0)) = 
l/\/2(|10) + 1 01) ) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 (dashed 
curve), 7 = 0.6 (dot dashed curve), 7 = 1.0 (dotted curve). Other parameters are (jlq = 2g, go = g, 
T = g. 




gt gt 

(a)C(i) 0>)Fd(t) 

FIG. 7: Time evolution for (a) Concurrence, (b) Fidelity from an initial two-qubit state of |^(0)) = 
l/v2(|10) + 1 01) ) at different values of temperature: T = 0.2g (solid curve), T = g (dashed curve), 
T = 5g (dot dashed curve). Other parameters are ^0 = 2<7, go = g, 7 = 0.2. 

and the fidelity gives the information of real state evolution. 

In Fig. [7] and Fig. [HI we plot the dynamics behavior of the concurrence and fidelity 
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gt gt 

(a)C(t) (b)Fd(t) 

FIG. 8: Time evolution for (a) Concurrence, (b) Fidelity from an initial two-qubit state of \ip(0)) = 
l/\/2(|ll) + |00)) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 (dashed 
curve), 7 = 0.6 (dot dashed curve), 7 = 1.0 (dotted curve). Other parameters are /xo = 2g, go = g, 
T = g. 




gt gt 

(a)C(t) (b)Fd(t) 



FIG. 9: Time evolution for (a) Concurrence, (b) Fidelity from an initial two-qubit state of \ip(0)) = 
l/\/2(|ll) + 1 00)) at different values of temperature: T = 0.2g (solid curve), T = g (dashed curve), 
T = 5g (dot dashed curve). Other parameters are ^0 = 2<?, go = g, 7 = 0.2. 

at different temperatures. When the temperature is as low as T = 0.2g, both cases of 
V>(0) = l/v^QlO) + |01> ) and ^(0) = l/y/2(\U) + |00>) display a periodical oscillation, 
the concurrence can always nearly restore its initial value. But for the former case, the 
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dynamics of fidelity is synchronous with that of concurrence; for the latter, the revival of 
the concurrence does not always mean that of the state. From the viewpoint of the definition 
of fidelity (I26p . it is partly due to the system part of the Hamiltonian (121). In the special case 
of the bell state l/v^GlO) + |01)), Pideai{t) is identical to the spins initial density matrix 
(It is an eigenstate of Hs), while for ip(0) = l/y^Gll) + 1 00) ) Pideaiif) is not in the same 
condition. It is further proved that the properties of the dynamics should be determined by 
the combination of concurrence and fidelity. 

V. CONCLUSION 

We have studied the dynamics evolution of two separated qubit spins in a bath consisted 
of infinite spins in a quantum anisotropic Heisenberg XY model. The bath can be treated 
effectively as a single pseudo-spin of N/2 spin degree. After the Holstein-Primakoff trans- 
formation, it will further be considered as a single-mode boson at the thermodynamic limit. 
The pair of qubits with no direct interaction served as an quantum open subsystem are 
initially prepared in a product state or a Bell state. Then the concurrence of the two qubits, 
the ^-component of one of the subsystem spins and the fidelity of the subsystem can be 
determined by a novel polynomial scheme during the temporal evolution. It is found that 

(i) larger anisotropic parameter 7 can help to maintain the initial state of the two qubits; 

(ii) the bath at higher temperature plays a strong destroy effect on the entanglement and 
coherence of the subsystem, so does the one with strong interaction g Q ; (iii) the dynam- 
ics of the subsystem is dependent on the initial state and in some special cases, only the 
concurrence is not sufficient to judge the revival of the subsystem. 
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